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CAWnjIAIICRS Of OHE HOOK'S HE&T HISTORY AT DIKFERJiNT COMCMlTiiATI-JNe; 

OP ll&SIOACTIVB SLBIl^TS TAKING ACCOUNT 0£P THE MATERIAL DlFFERiUTIATICJS 


NITH UBI/TSIG 



0 . 1 o Oma'tak^a , X a « 1 . A1 * ber , I . LoHy aaants e va 
I «. Introduction 

Papers by GcUrey,G.MacDonald,2*A,I.ubii!iove,BoXuiievin^S. V oMaeva^ 
PcPriAer^H^Heynolda and AoSuojtters, J®Iriyama euid J«Schimazu,D. Anders cm 

andRcPhinerey jT®Hanks and D®Anderson.,R,MacGonnelsH«Tokso2 and SeSol omon 
et al B'"I7f265 29*'3ll are devoted to tb.e investigation of 

tlie Moon*s h,eat history based on the solution pf the heat conductiv:i t.> 
equation with the given values of the Hoon*s material parameters^ the 
initial and boundary conditions and with the heat sources-^ radioactive 
elements with the given concentrationsc Our papers [20,2l]are devoted 
to the same problem® Up to the homogeneous and laminated 

Moon models were considered .separately® Then in 1966-6? we and 
pective of us P®Prik8r,R#Reynolds and A. Somers, took jnfco account 
the calculation the differentiation of the Moon* s material and convt!c-» 
tion of the radioactive elements to the surf ac e with melting. The con- 
sideration of continuous differentiation in the calculation peraitj to 
investigate the Moon*s heat history in more detail. In I9?i the dif- 


ferentiation was also taken into account in the paper by S.V®Maeva [13] 
and in 1973 in the paper by HaToksoz and S.Solomon [26.1 


Thusj presently j the calculation scheme of the Moon's heat his-- 
tory is sufficiently full developed by the computer and much new re- 
sults are obtained^ fet this problem is of Interest so far and this 
interest stiii increases. This associated w.th the large uncertaxaty 
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in the !i4oon*s parameters values enter jLog ia. the equatioa of the heet 
conductivity • This reTers also to the uncertainty in the given vaXuee 
of the density* the heat capacity and the heat conductivity and their 
dependence on the Moon's material state (for eacaaplet idMn melting) 
and to the uncertainty In the given curve of melting to . the initi<> 
al conditions (in particular* to the Initial temperature) and to the 
Moon's age related to Its origin* However* the determiziatlve in the 
Moon's heat hlstdry is the accepted concentration of radioaotlfe eje- 
ments in the Moon's material and their redistribution In the Moem 
interior in time* 

In the above paper except fo [20*21] ttio concentration of radio* 
active elements Is tahen according to their contents In chondrlte me* 
teorites (for example* Urey azul MacDonald}* or in '^e Sarth (Levin* 
Maeva)*or In the "Earth mixture". The values of the heat flux through 
the Moon's surface were obtained to he of ^e order of (0«3*0.6) 10*^ 
cal.cm“*^sec"^. 

At the same time the results of radio astronomical InveatigatlonB 
by V.S.Troitaky and V.D.£roti^ov (27f2Si give the fl^ value equal 
to 0.35-0.95 10“^ cal cm^^sec*"^ which was considered for a long time 
to be rather doubtful but was confirmed by the resent intaatlgatlons 
of the lunar material frcm Appolo*I5 t?!* is a result It is necessa^ 
to calculate the Moon's heat history at different variant of the 
radioactive elements contents which may account for so large heat 
flux. 

The first attempt of this calculation was made in our papers in 
1966*1967 [20*2x1 . Ve use four variants of radioactive elaments con^ 
tents two of which corresponds to concentrations accepted by V.Xu. 
LeV’in*S.V.Maeva [9I (variant "C mean contents in the Sarth) and 
G.Urey and G.Mac Donald, and the rest two give relatively large con- 
centrations with the flux value close to the radio a8troncmiical|L data* 



The assumptioa was \ised on the exponential differentiation (when 
meltlzis) of radloactlYe eleaezits* Since radioactive elements with 
large concentration were close to the surface of the Moon (in I0-2Qkm) 
they weakly influence dn the temperature of the main liu^r mass which 
quickly cooled* As a result the thickness of the hard cover in this 
variants was equal to 600«>700km, l.e* larger than at the mean (by B* 
Yu*Xi6Vin) and at the small (by Urey « MacDonald) concentratlona where 
the thickness of the cover was 230km and 400 km« 

We use in C20»2llthe exponential model of radioactive elements 
covection to the surface at differentiation corresponds to a certain 
degree to the so»colled *' sudden diff erentiatlon”. This model permits 
to find out the main physical peculiarities of the heat history rela- 
ted to differentiation yet It limits the number of variated parame- 
ters* In this case laminated differentiation proposed by Friker,et al 
5 ] is proven to be evidently close to the real processes- That is why 
we used in the present paper the laminated differentiation of radio- 
active elements when melting with different portions of elements 
(from 100% up to 20%) convacted from layer to layer. The variants of 
concentratic^ of radioactive eluents were taken the same as in 12.0^ 

21.1 

It should be noted that we are not specialists in geophyslk and 
we have no rights to discuss questioxis of the Moon origin and its 
chemical composition* For this reason for the calculation we take 
parameters either from the expected values of the heat fluxes (at 
the given concentration of radioactive elements) or from the fundamen- 
tal papers by other authors. Thus the initial temperatures (273^^ 
aad 3QO°K) —re taken on the ba^ia of papers 6,6] , the beat capaci- 
ty ” tbe basis of tbe paper by S«V»Maeva wbere a sbarp* 

change of the beat capacity at the time of melting is considered. 
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Remarics bj Idib imova (X2] are taken into account on tlje possibility 
of the flux increasin^^ at the reasonable concentrations of radioactive 
eiemnts considering the contribution from exciton coBpoiieat of the 
heat conductivity and froin the amplification of the efficiency of tte 
radiatioaal component. 

Consideration carried out Involves significant idealization ref- 
fered first of all to the process of l^ination* That is why the re 
suits obtaines are rightful only in the limits of our mathematical 
model. 


The statement of the problem 


As a rule the Moon is assumed to be a sphere with the radius 
^ =1735 km, with homogeneous density p =3»3^/cm^, with the heat 
capacity C and the heat conductivity K, defined according to 
Lub imova hs] (heat conductivity ) and Maeva (heat capacity) 

by the following relations: 

The best capacity 


C(T, f^(T)] 




* • ttfcJTi 




3— 


L =100cal/gm Is the heat of meltings s200 £; 

= in cm, T = in K. 

+ 500 (I - ) is tlie begimi&g of tUe sate- 

rial melting, of the silica (Uray [29l ). 

We consider in (I ) gradual absorption of meltixsg heat In the 
interval of melting temperature aT (the effective heat capacity in 
the interval aT first increases linear and then decreases) • 
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The heat coxiductivlty 


K= i(T) + 




3€ 




( 2 ) 


In (2} the grld( phonon) , radifttlonal (photon) and ezciton components 

of the heat conductivity Cxz] are taken into account. 

/ 

Phonon heat conductivity outside the interval of melting chaiages 
Inversly proportional to the temperature and during melting It shaply 


increases'*'^ and remains constant: 




9 CaB. 


^(T) = 0 . 0 ^ 


cae 


C/rt. 




cm., y^at 
s con^ , 




Photon and exciton components of the heat conductivity are defined 
by the well known relations where Xs the refractive index, 

0 =4.27 10”^ cal/cm grad year is Stephan-Boltzmann constant, 
^=25 cm"*^ is the absorption coefficient of the material, f =5il7^ 
I0“^^cal is the energy of exciton activatioi'^\ 

In the assumption of isotropic distribution of the heat sources 
and when the beginning of coordinate is in the centre of the sphere 
the heat conductivity equation is written in the forms 


/'*’ '3'* 




(5) 


where H is the heat sources «> long lived radioactive elements * 
uranium, thorium, kalium-40. 


^See the paper by Tikhonova^ Lubimova, Vlasova 

“*"*‘^The values are taken for the pyrolit (the mixture of basalt and 
peridotit) ^Moon's composition similar to the Earth mantle - see 
Ringwood [25] • 



^ ^ Q ' . ( 4 ) 

J " 

where the fuactioa def taes the distributioa of radioactive ele- 

ments depending on the radius a£ different time. Indexes j =1,2,3 
refer to uranium, thorium and icai ium-40, respect ively,J^= I *54 IO"J§ar*^ 

A =5 I0~^^ year-^ , yj, =5.7 I0”^° year"^, the heat release of one 
gramm of the radioactive element is 0*805 cal 

oCj^=OeI93f 0(35=0*224 and values of - the concentration at the 

beginning of consideration - are defined based on the expected presen- 
tly concentration of over the formula Hj® -r 

where t^=4*7 10^ years is the Moon age. We consider four variants 
of values of radioactive elements concentration given in table X. 
Radioactive elements concentrations are chosen based on the stationary 
heat fluxes equal to 1. 35 and O.91 cal cm*“^sec~^ for variants I and 
II and corresponding to radio astronomical data by Troitsky and 
Krotikov f27j and to the Appolo-I5 ffj , and 0.61 and 0.236 cal.cm**^ 
sec for variants II and III coinsided with one of the variants by 
Levin and Maeva f9l (II) and close to concentrations accepted by 
Urey and MacDonald (III). The relation of thoritim concentration to 
uranium one is taken to be equal to 4 (see [*93 )> uranium concentra- 
tion is taken in variants 1,1'*’^, II based on the pyroiit Moon compo- 
sition and chondrite composition in variant III. This assumption on 
the pyroiit Moon composition is founded on the similarity of the 
Moon composition and the Earth mantle noted by Andez*son fl3 (for 
example, the mean atomic weight of the Earth mantle M =22.4 which 
refers to ejected rocks, coaly, ordinary and enatative chondrices. The 
same atomic weight has the Moon* s . material 22.0) axid also on the 
data by Sorveyer - ? and on the paper by MacCrea fl^fi^eyaolds, Summers 
f243 et al, i?rom all these data it follows that the Moon consists of 
sil:!.cate (for example, basalt) similar to those at the Earth 


c 


- 7 - 


xoantle'*'^ )* The table gives uranium concentrations and the relation 
of K-4-0 and uranium concentrations for different rocks* 


Hocks: 

Grqnite: 

Basalt: 

r* 

Ak logit e: 

Peridot it I 

Dunltt 

Chondrite: 

U-IO"® 

400 

80 

4.3 

0.6 

O.I 

I.I. 

K40/U 

1.05 

1. 13 

1.44 

1.2 

1-2 

9 


From this table and the t^lven stationary fluxes «e choose the above 
variants and II with the content of basalt equal to 14%, 10% • 

and 3%f peridot it ~75%,77% and 82% and Iron ~ 13% (Reynolds, Summers 
f24] ) with impoverished content of radioactive elements (Iron, has 
U=0.I-I0~®g/g, Th=0.4*I0"®, =0) and the ralatloa k'‘®/0sI.I3 

and 1.6 (see Levin [lO]and Lubimova [X2] )• We take also variant III 
- chondrite model (Urey and MacDonald), 

rJote that variant II corresponds to variant ’*€*• of the mean con- 
tent of radioactive elements in the Earth (see Staxkova C333and. I*u- 
bimova fl 23 

As the initial condition we take two variants of the "Cold” 
and the ”hot" Moon*s model at the moment of its formation according 
to initial temperatures 2?3^ and 900^. These are the same initial 

condition taken by Urey, MacDonald, Iri 3 ama,Schima 2 u,Priker'*"*‘'*'^ based 

■ ' ^ ' 

^^Prom this it is not necessary followed that the hypothesis is true 
on the Moon formation by the way of tearing off the Earth since even 
interstellar material contains grains of high melting silicates 
(Mac Crea (193 )• 

^"^^Runniiig ahead we note that just this variant gives the flux corres- 
ponding to Appolo-i5 data and the radio astronomical one8,l*e» eviden- 
tly. the terrestrial composition of the Moon Is rather probable. 

■♦■++) r 1 

Maeva took the parabolic tenperature distribution with the 

value in the Moon’s centre equal to 500^. 



« o 


on. the di_fferent n;ypoth,esis of the Moon origin and the tine of Its 
accumulation and the initial heating we do not give dona Ide-at Ion. 
The bjundary condition ( ) is defined by relation (5) 
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(5) 


wherefis 0.9 I0**^cal cat ^sec ^ conslderea the Moon heating by the 
solar radiation. This condition leads practlcalljr at the calcu\a*^ion 

to the sane results that the more simple condition T » 273^ sconst 
does. 

The fflechanlam of differentiation 

The mecoanism of diff erentiaticm is defined by the function of 
distribution of radioactive elements F(^»i) which is given by the 
following manner. At first, we consider a homogeneous model of dis- 
tribution of radioactive elements, Ue. P(r^i)\x (analogous to 
f20,2I ] ) is the first stage of calculation. Then, after reaching at 
the moment , at the point , the temperature corresponding to 

the beginning of melting T^^^^ the distribution of tempera- 
ture obtained to the time of (where is the time of 

melting — see below} is taken for the new initial distribution, 
the distribution of radioactive elements is assumed to be laminated» 
namely, radioactive elements from the layer transfer to 

the adjacent upper la*^er ( i ) and the calculation is carried 
out with the poored layer ( However, for the time At de- 

fined by the necessity of the layer heatingby by aT ss2D0\ more 
deep layers have reached the temperature of meltlzsg (since melt 

increases with the depth, x.e. with the decrease of ). JRadioactlve 
elements from these layers go to the layer , i.e. their 

content in this layer is equalized and so on, until the lower layers 
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stop melting. Such laminated differentiation (similar to there accep- 
ted by F rick er, Roy no Ids, Summers) is defined by five variants differed 
by the value fl waich is the portion of radioactive elements con- 
vection from the layer* We accept fl =I (the complete convection), 

0*a (80% of radioactive elements convected from the layer) ,0. 5,0.4, 

0.2, In this way after difi'erentiation in each of the melted layer 
there remains from 0 up to 80^5 from the content of radioactivj elements 
in nondiff erential material (note that in the latest paper by Maeva 
[i5] the residue was from 2% up to 20%,i.e* 0.8). 

The present paper does not consider the possibility of convection 
of the melted material to the external hard layer as it was done jn the 
paper by Maeva [15J since the concentration of radioactive elements 
in the thin surface-attached layer leading, in reality, to the lore 
rapid cooling of the Moon and to the larger thickness of the hard 
cover is obtained in the model of exponential diff erentiation which 
is considered in detail in [20, 2l] , 

3* The method of calculation 


For the numerical solution of the heat conductivity equation 
we use the method of line coordinates where the section • 

(where Ha is the radius of the Moon) is divided into 100 l§yers 
and only arbitraries over the special coordinate are replayed by the 
difference ratio. Thus the solution of the boundary problem is reduced 
to the solution of Cauchy problem for the system of ordinary differen- 
tial equations. It should be noted that the approximation error the 
larger the closer is the layer to the surface. That is why nonuaiform 
grid is used in the calculations the layer thickness close to the 



-- 10 - 


surface ( ^ 5km) la in several times smaller In comparison with the 
depth layer ( ^ 18 fcn). The system obtained is solved by the appro- 
ximation mefchod providing the accuracy of temperature calculation 
(at each step over the time) of the order of 10*^ « All calculations 
are made by the computer B£SM -6 • 

4. Results of calculation 

Results of calculation are given in ?lgs I-X4 arl in Tables 2 
and 3* 

The process of melting and differentiation 

Figs 1-4 shows the time temperature dependence in different 
Moon layers (along the radius beginning from the centre 0 ) for all 
variants at fl si • Here also temperatures are plotted of the beginning 
of melting corresponding to each layer and the time dependence of 
the beginning of melting in different layers* From Figs 1-4 it is 
claer that melting iu all variants begins in the layer near 

ISOO km from the centre (l*e* at the depth ^ 250 km under the sur- 
face) at the time period from 0*4 up to 1*5 milliard years from the 
beginning of the Moon history) and rapidly (for O.i -0*5 milliard 
of years) the region of melting propagates to the centre and to 
the Moon*s surface* At the same layer where melting begins for the 
first time, in 0.2 -0*5 milliard of years the differentiation of 


‘'‘^Altogether there are 100 layers. 



II - 


the material begins for the first time. The region of melting at the 
moment of the beginning of differentiation "^iixit.dif. i^eaches 1600km 
and at the same time a sharp temperature increase occurs in the 
layers close to the surface. As a result the temperature in these 
layers is close to or exceeds melt. • It is clear since 

melted radioactive elements convected to these layers from the more 
deep layers'* The deep lasers (the depth more 500«-600 km,i*e. 
)^s:1100-1200 km) to the moment of differentiation turn to be heated 
higher than -iai-t.jBelt slowly cooling due to the convection 

of radioactive elemants to the upper layers'**^. From Pigs 5-7 it is 
seen that in this way almost simrltaneous heating of the Uoon mate- 
rial occurs up to the depth of 200-250 km \mder the surface. Layers 
which are closer to the surface heating slower and the depth of these 

layers increases as time goes by. Numerical data are given in 

^ •* 

Table 2 and 3» 

The beginnlriK and the duration of the maximum melting of the interior . 

Figs 8-XX give time dependence of the radius of melting. From 
these Figures it is clear that in all variants the period of the 
maximum melting is observed continued 1-2 milliard of years (depen- 
ding on the variant - see table 2 and 3 ) beginning In the inter- 
val 0.7-2. 3 milliard of years from the beginning of the Moon history. 
It is naturally that for large concentrations and **hot” models the 
maximum melting begins earlier than for the ooncentrations of small 
and ”cold” models. So for the. variants I and I"*^ the maximum melting 

'‘^Differentiation occurs in all variants except for III (cold model) 
where the energy of radioactive elements is not sufficient for th€ 
beginning of differentiation. 



- 12 - 


bee ins approximately 5«5 - ^ milliard of years ego aad coatinues 
aurlng 1.5 “ I milliard of year *5 and in the case rf small (variant III) 
end mean concentrations and the large convection (fl^O^S) of radio- 
active elements (variantll) the maxinwm melting begins 2*5 - 3 aii- 
liard of years ago and continues I»5 - ^ milliard of yearst x.e. it 
begins later and continues longer as should be 'ucpected'*’^ The pecu- 
liar case is variant II (’’cold’' model) at the siauH convection ra- 
dioactive elements ( Jl^O.6) vhere the period of maximum melti ng be- 
ginning 2.5 milliard of years ago and continues at the present time. 

Here the specific combination of conditions Is realised^ The concen- 
tration of radioactive elements is not so large end the process of 
differentiation is not so effective that a fast convection of radiD- 
active elements to the surface occurs when their influence on the 
further hit^tory turns to be small but not sufficiently small for the 
natural cooling of the Moon material* 

Thus, in all variants except for the mentioned II ( fi^O.6) the 
cooling of the Moon takes place presently* 

The thickness of the bard cover 

As it *is seen from Figs 8— II and Tables 2 and 3 the t hi c k ness of the 
hard cover at the period of maximum melting ranges in variants Ifl^ > 
II, III 15-20 km, 25-30 km, 35-^5 km and 100-200 km , respectively, l.e. 
up to I - 2.5 milliard of years ago the Moon was almost completely 
melted (in variants II ^ith rv^0«6 the Moon even now should be mel- 
ted up to the depth 45 km under the surface), then cooling beg^ 
and presently the thickness of the bard cover in all variants is 
150-250 km (except for the mentioned "peculiar caae"}* The same is 

'*’‘^The smaller concebtration is the more quite pve the process passed. 



followed from Figs 13—12 ■,7iiere th3 distribution of the temperature 
in the Moon at the present time is given* 

The heat flux through the Moon surface 

Fig 13 gives time dependence of the heat flux through the lunar 
surface* From these ctirves it is clear that at the period of the 
maxi mu m melting of the Moon flux as it could be expected is also 
aaximm one exceeding 2-3 times the stationary value the given 

concentration of radioactive elements* Now at all variancs (except 
for II with the flux decrease is observed and its approxima- 
tion to the stationaiy value >!^©t even now the flux exceeds 

nearly I.S times • The flux value in variants II 

(h.9^0.8) and III is equal to 1*8 *10“^ cal cm“^sec~^, 1*27*10"^, 
and 0*9 - Oii95*IO cal* cm^^sec*"*^, respectively, in variant III 
(’’cold” and ”hot** models) - 0..28 and 0*33 I0“^cal cm^^sec"^ and In 
variant II (H. ^ 0.6) is a. rather large flux equal to 1*6 I0"^cal.cm“^ 
sec that associated with the melting in this variant extended up 
to the depth of 43 1cm* 

Thus fluxes close to the measiired ones by the radio astronomical 
methods and obtained by Appolo-I3 are realized in variants I"*”** and 
II (R^ 0»8) ,i*e« at the terrestrial or “mantle” - pyrolit content 
of the Moon material* 

5e Conclusion 

Calculations are given where different cone ehtrat ions of radio- 
active elements are accepted which refered both to the chondrite 
composition (III) and to the terrestrial one (II) an d to the Barth 
mantle -pyrolite composition (I and 1-^+). The material differentia- 
tion and the convection of radioactive elements to one surface due to 
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in jlting are takea iato account at the process of calculation over 
the scheme close to the scheme by Fricker^ReyaoZds and Summers, 
Calculatioiid show that the Moon*s interior has been heated up to the 
melting during the first 0.7-2«5 miJliard of years after its formation. 
The maximum melting of the Moon covering practically the entire 
Moon up to the depth 15-45 km from its suifaco began 5*5-4 milliard 
of years ago in the assumption on the pyrolite composition of the 
Moon or 2,5“5 milliard of yejrs ago in the assumption on the terres- 
trial or chondrite composition of the Moon, This maximum melting 
continues during 1-2 milliard of years and at present the Moon is 
cooling (there is a scarcely probable variant of weak differentiation 
at the "terrestrial*’ composition of the Moon when melting is conti- 
nuoing even now). The present depth of the hard co\er is T50-250 km 
and the value of the heat flux approximately 1*5 times exceeds the 
stationary one, and is equal for variants I, I'*'**’, II and III to I# 8, 

1. 27 » 0,9-0, 95, and 0,28-0.55 10"^ cal cm“^sec~^, respectively. 

From Tables 2 and 5 it lollows that in the most probable vari- 
ants with pyrolite-^mantle” (I^'*’) and terrestrial (II) composition 
of the initially "hot" Moon the maximum melting of the Moon began 
5, 5-3, 3 milliards of years ago thet corresponds to the age of the 
basalt and mask ones (3* 13-5. 71 milliard of years). This testify to 
the vulcanic activity and was completed approximately 2,7-3 milliard 
of years ago that agrees with the absence of the material samples 
at the Moon yunger than 3 milliard of years (Papanastassiou,Wasser- 
burg C22,23,32j ) when the activity stops« So, limitations stated 
resent ly by Toksoz and Solomon f26) are satisfied| yet the thick- 
ness of the hard cover obtained by us is essentially small (In £26j 
the thickness of the cover is 600 km) , 

It should be noted that the most close to reality is the vari- 
ant of the "terrestrial** composition of the Moon (composition **C** 
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according to Stark ova) giving inspite of jnoderate concentration 
a sufficiently large heat flux'*'^ corresponding to the radio astrono- 
mical measurements and Appolo-I5 data, i*e« to obtain the large 
flux observed it is not nocessary large concentrations of radio ac- 
tive elements to be present. 

The authors are grateful to Prof .V.S.Troitsky and i'rof.N.L'c 
Tseitlin for the attention to the present paper and for useful dis- 
cussions and to M.V.Ostrovskaya for the help in the construction 
of numerous graphics* 




This possibility was noted by Lubimova In fisj 
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